% Correr desde octave 
% lab5ej2(N)
% Donde N es el numero de intervalos: 
% 4,10,20,etc
function lab5ej2(N)
		
    if nargin == 0
        printf('\nError, ingrese "N" numero de intervalos: lab5ej2(N) \n');
        exit;
    endif
    
	a = 0;
    b = 1;
    fun = 'f';
    Y = quad(fun,a,b);    
    
    printf('\nCalculo de area para funcion exp(-x) entre 0 y 1 con %i intervalos\n',N);
    
    printf('\n_________ Trapecio compuesta __________\n');
    
    regla  = 'trapecio';
    
    S = intenumcomp(fun,a,b,N,regla);
    printf('\n intenumcom: 	\t%f', S);   
        
    printf('\n integral exacta: \t%f', Y);
    
    printf('\n error absoluto:  \t%f \n' , abs(S-Y))
    
    
	printf('\n_________  punto medio compuesta _______\n');
    
    regla = 'pm';
    
    S = intenumcomp(fun,a,b,N,regla);
	printf('\n intenumcom: 	\t%f', S);   
    
    printf('\n integral exacta: \t%f', Y);
    
    printf('\n error absoluto:  \t%f \n' , abs(S-Y))
    
    printf('\n_________  simpson compuesta __________\n');
    
    regla = 'simpson';
    
    S = intenumcomp(fun,a,b,N,regla);
    printf('\n intenumcom: 	\t%f', S);   
    
    printf('\n integral exacta: \t%f', Y);
    
    printf('\n error absoluto:  \t%f \n' , abs(S-Y))
        
    printf('\n');
    
endfunction


function y = f(x)
    y = exp(-x);
endfunction
